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Inelastic light (Raman) scattering has been used to study electronic excitations and phonon anomalies in de- 
twinned, slightly overdoped YBa2Cu306.gs and moderately overdoped Y0.85Ca0.15Ba2Cu3O6.95 single crys- 
tals. In both samples modifications of the electronic pair-breaking peaks when interchanging the a- and 6-axis 
were observed. The lineshapes of several phonon modes involving plane and apical oxygen vibrations exhibit 
pronounced anisotropies with respect to the incident and scattered light field configurations. Based on a theoret- 
ical model that takes both electronic and phononic contributions to the Raman spectra into account, we attribute 
the anisotropy of the superconductivity-induced changes in the phonon lineshapes to a small s-wave admixture 
to the d x 2_ y 2 pair wave-function. Our theory allows us to disentangle the electronic Raman signal from the 
phononic part and to identify corresponding interference terms. We argue that the Raman spectra are consistent 
with an s-wave admixture with an upper limit of 20 percent. 



PACS numbers: 74.20.Rp, 74.25.Gz, 74.25.Kc, 74.72.Bk 
I. INTRODUCTION 

The symmetry of the pair wavefunction provides impor- 
tant clues to the mechanism of high temperature supercon- 
ductivity. Phase sensitive experiments have shown that the 
dominant component of the order parameter exhibits d x 2_ y 2 
symmetry^ In cuprates of the YBa2Cu3C>7-type with or- 
thorhombic crystal symmetry, an s-wave admixture result- 
ing in an anisotropy of the superconducting (SC) energy gap, 
2A, between the a- and &-axes is expected. A more gen- 
eral trend of increasing s-wave admixture with increasing 
doping level has also been reported] 3 ' 4 ' 5 ' 6 In YBa2Cu307_5 
( YBC07_a), such an admixture has been confirmed by angle- 
resolved photoemission spectroscopy (ARPES)^ and flux 
measurements on Josephson junctions^ However, the inter- 
pretation of the ARPES data is complicated because the sur- 
face properties of YBC07_,5 may differ from those of the 
bulki^ The phase-sensitive measurements, on the other hand, 
have yielded only lower bounds (9%) to the s-wave admix- 
ture. In addition, phenomenological calculations based on 
the two-dimensional Hubbard mode l 10 ' 11 have demonstrated 
that the in-plane anisotropy of the spin fluctuation spec- 
trum determined by inelastic neutron scattering on detwinned 
YBCO 12,13 is consistent with a small isotropic s-wave ad- 
mixture to the d x 2_ y 2 -pairing symmetry. However, the same 
calculations show that the dispersion and spectral weight of 
the spin fluctuations are also influenced by anisotropic hop- 
ping parameters that are difficult to determine independently. 
Moreover, different interpretations of the neutron data have 
also been propose d 1415 . Further experiments are therefore re- 
quired to conclusively establish the magnitude of the s-wave 
contribution to the SC gap in YBCC>7_,5. 

Inelastic light (Raman) scattering is a powerful probe of 
electronic and lattice vibrational excitations in high-T c super- 
conductors. 16,1 7 Both electronic and phononic Raman scat- 
tering have been applied to study the superconducting gap 
anisotropy in YBC07_^. The anisotropy can be inferred 



from the energies of gap features of the electronic contin- 
uum in various polarization geometries* 6 *^ but the accu- 
racy of this method is limited by overlap with Raman-active 
phonons. Possible manifestations of the gap anisotropy and 
an anisotropic electron-phonon coupling in the lineshape of 
a Raman-active out-of-plane vibration of the in-plane oxygen 
atoms of Big symmetry have also been investigated^ How- 
ever, the analysis of these data has been contested- - and a 
quantitative estimate of the gap anisotropy has not been ex- 
tracted from them. 

The present work was in part motivated by a recent theo- 
retical study that yielded quantitative predictions for the elec- 
tronic Raman continua in a (d x 2_ y 2 + s)-wave superconduc- 
tor^ In order to enable a detailed comparison with these pre- 
dictions, we have performed Raman scattering measurements 
of the temperature dependence of the electronic continua and 
phonon modes in twin-free, slightly overdoped YBCC%.95 and 
moderately overdoped Y0.85Ca0.15Ba2Cu3O6.95 (henceforth 
YBCC>6.95:Ca) single crystals. We found that interference 
between the electronic and phononic scattering channels im- 
posed severe limitations on our capability to extract infor- 
mation about the gap anisotropy from the electronic contin- 
uum alone, in accordance with previous worki 6 ^ We have 
therefore employed a phenomenological model that treats the 
electronic and phononic contributions to the Raman signal 
on equal footing. This formalism yields predictions for the 
energy, intrinsic linewidth, and Fano parameters of phonons 
coupled to the electronic continuum as a function of temper- 
ature. From a comparison to the experimentally determined 
lineshapes of two phonon modes, the Bi g vibration of the in- 
plane oxygen at 340 cm -1 as well as an apical-oxygen vibra- 
tion at 501 cm -1 , we extract an upper bound of 20% on the s- 
wave admixture to the SC energy gap. This is consistent with a 
10% to 15% admixture of s-wave to the d x 2_ y 2 pair-breaking 
peak observed in the electronic continuum in the supercon- 
ducting state, and it coincides with the lower bound extracted 
from the phase-sensitive measurements^. Similar effects have 
been recently reported for an organic superconductor with the 
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same point group as YBa2Cu3C<7_,5 (£*2/i)- In this case, how- 
ever, the analysis led to the conclusion that only an isotropic 
s-wave state is present.— 

This paper is organized as follows. In Sec. [ID we describe 
experimental details and discuss the raw Raman spectra. An 
introduction to our theoretical model is given in Sec. [HI] In 
particular, we focus on a simultaneous description of both, the 
electronic Raman response and the phonon anomalies above 
and below T c . Based on this description we are able to an- 
alyze the superconductivity-induced changes in the phonon 
lineshape in Sec. [IV] A summary of our results and a com- 
parison to prior work is contained in Sec. [V] 



II. EXPERIMENTAL DETAILS 

High quality single crystals of YBCO6.95 and 
YBCOe.95:Ca were grown by the top-seeded solution 
growth method.? 3 - The orientations of the crystallographic 
axes were determined by Laue x-ray diffraction and polarized 
light microscopy. The samples were then cut into rectangular 
shapes of typical size 3x3x1 mm 3 , annealed at 520°C in a 
flow of oxygen gas for 150 hours, and quenched in liquid 
nitrogen in order to avoid further oxygen diffusion. The mag- 
netization curves of the YBCO6.95 (YBCC>6.95:Ca) crystal 
show an onset of T c = 92 K (75 K). A transition width AT C 
of less than 3 K (6 K) indicates good homogeneity of the sam- 
ples. Employing Tallon's phenomenological expression,? 4 
the hole doping level of our samples is estimated as 0.17 
(0.21) per planar Cu ion, i.e., they are slightly (moderately) 
overdoped. The crystals were detwinned individually using 
the procedure described in Refs. |25| andU^. Raman spectra 
measured at different spots of the afo-surface of each sample 
look basically the same, thus confirming the homogeneity of 
our single crystals. 

The Raman scattering experiments were performed in 
backscattering geometry using a Dilor XY-triple grating Ra- 
man spectrometer with a CCD camera as detector. For exci- 
tation, the 514.5 nm line of an Ar + /Kr + mixed gas laser was 
used. The resolution of our spectrometers for this experiment 
was about 3 cm" 1 . In order to avoid heating of the sample, 
the power of the incident laser was kept below 10 mW at the 
sample surface with a laser spot of 100 /m in diameter. The 
direction of the incident laser light was always parallel to the 
crystallographic c-axis. For the selection rule analysis, a po- 
larizer and an analyzer were placed into the light path before 
and after the sample. 

YBCO6.95 and YBCOe.95:Ca are orthorhombic with the 
space group Pmmm (-D271 symmetry). In polarized Raman 
scattering experiments, excitations of A g , B\ g , A g +B\ g , and 
A g symmetries are accessible by using xx (or yy), xy (or yx), 
x'x' (or y'y'), and x'y' (or y'x') polarizations of the incident 
and scattered light fields, respectively»ii2i In this notation x 
and y correspond to the direction of the electric field of light 
along the a- and 6-axes, whereas x' and y' correspond to the 
diagonal directions, i.e. x' ~ x + y and y' ~ x — y, re- 
spectively. As most previous work on YBCO has been per- 
formed on twinned specimens, we follow previous publica- 
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FIG. 1: Raman spectra of a detwinned slightly overdoped 
YBa2Cua06.95 single crystal in xx and yy polarization taken with 
an Ar + laser line (A=514.5 nm). The mode assignment corresponds 
to Refs. [§2.. 33 3435.36. The spectra were shifted by a constant off- 
set with respect to each other. The intensity scales in left and right 
panels are the same 



tions and use the tetragonal notation of the polarization sym- 
metries throughout this paper. In this notation, the structure 
is described in terms of the closely related tetragonal point 
group £>4h, and the excitations of A g +B\ g symmetry are ac- 
cessible for either xx or yy polarization, B2 g for either xy or 
yx polarization, A g +B 2g for either x'x' or y'y' polarization, 
and Bi g for either x'y' or y'x' polarization of the incident and 
scattered electric field vectors. 

Fi gure [T] show Raman spectra of detwinned YBCO6.95 
crystals measured at various temperatures. As in other high- 
T c superconductors ^ 16 i 28 - 29 i 30 i 3 1 me normal-state spectra ex- 
hibit a flat electronic continuum with superimposed phonons. 
Below T c , a significant fraction of the electronic spectral 
weight is transferred to higher energies, resulting in a broad 
pair-breaking peak in the electronic continuum. In addition, 
the phonon lineshapes reveal characteristic changes. In Bi g 
polarization this spectral-weight redistribution is most pro- 
nounced because electronic Raman scattering in this geometry 
is sensitive to the maximum of the SC gap along the anti-nodal 
direction of the two-dimensional SC gap. We therefore focus 
in this work on this geometry, along with the xx and yy po- 
larization channels that provide direct information about the 
in-plane anisotropy of the electronic response. 

Figure Q] gives an overview of the phonon modes of 
YBCO6.95 in the xx and yy polarization geometries. As ex- 
pected based on a group-theoretical analysis , 17 - 32 i 33 i 34 - 35 i 36 fi ve 
phonons appear in both polarization channels. The lowest- 
frequency phonon at 113 cm -1 (Ai g ) originates predomi- 
nantly from vibrations of the barium atoms. The next low- 
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est frequency phonon at 148 cm -1 (Ai g ) corresponds mainly 
to the vibrations of copper (Cu2 ions). The vibrations of the 
apical oxygen ions (04) appears at 501 cm" 1 (Ai g ). The 
two modes at 340 cm" 1 (£>i g ) and 446 cm" 1 (Ai g ) originate 
from out-of-phase and in-phase vibrations of the planar oxy- 
gen ions (02 and 03), respectively. 32 Two additional defect- 
induced modes appear in the yy symmetry at 232 cm" 1 and 
at 579 cm" 1 . They originate from vibrations of the copper 
(Cul) and the oxygen (Ol) ions, respectively, in the Cu-O 
chains, which are aligned along the crystallographic 6-axisi 3 ^ 
In YBCO crystals with fully oxygenated chains, these two 
modes are Raman forbidden but infrared-allowed (Bi„ sym- 
metry)] 37 ' 38 ' 39 They become Raman-active due to breaking of 
the translational symmetry by defects, i.e., unoccupied oxy- 
gen positions in the Cul-Ol chains . 33 ' 34 ' 35 The absence of the 
strong mode at 579 cm" 1 in the data with xx polarization con- 
firms the high detwinning ratio of our crystal (~ 95%). 



III. ANALYSIS OF RAMAN SPECTRA 

The measured Raman intensity I a (uj) in a given polariza- 
tion channel a is related to the imaginary part of the response 
function Im% ff (w) via I&(uj) = 4[l + n(w)]Imx ff (w), 
where n(uj) denotes the Bose distribution and A is a coupling 
constant. The Raman response in cuprate superconductors in 
the optimally and overdoped regimes consists of electronic 
(intraband) and phononic excitations. In order to disentan- 
gle these two contributions it is necessary to employ a proper 
fitting procedure. A phonon interacting with intraband ex- 
citations acquires a renormalized self-energy and exhibits an 
asymmetric, Fano-type Raman lineshape. To describe such an 
electron-phonon coupled Raman spectrum the following for- 
mula can be employe d 40 ' 41 ' 42,43 ' 44 ' 45 



ImXff(w) = p a (uj) + 



9l 



T(uj) [l + e»] 



(1) 



x {S 2 (uj) - 2e(uj)S(uj) P<y (uj) - p 2 a (uj)} , 



where e(uj) = (uj 2 - n 2 )/2uj T(uj) and S(uj) = S + R<t(uj). 
The renormalized phonon frequency and the renormalized 



phonon line width are given by fr 



2uj glR a {uj), and 



T(uj) = To + g 2 p a (uj), respectively, with the intrinsic phonon 
frequency uiq and the bare phonon line width Tq. The param- 
eter So can be expressed in terms of the electron-phonon cou- 
pling g a , the Raman phonon matrix element T a , and the Ra- 
man electronic matrix element that is, So = T a j ( r y a .g a ). 
In our model, g a , T a , and 7 CT are assumed to be real and there- 
fore So is real. The Kramers-Kronig related functions R a (uS) 
and Po-{ui) denote real and imaginary parts of the electronic 
response function, respectively, Xa(^) = R<t(w) + ip a (uj). 
While the first term in Eq. ([1), p<j(uj), describes the "bare" 
electronic Raman response, the second term represents the 
phononic contribution and its coupling to the electronic back- 
ground. We will use expression (Q]i to describe the coupling 
of the electronic background to the phonons whose lineshapes 
exhibit the clearest manifestations of the electron-phonon in- 
teraction (e.g., the in-plane B\ g phonon). The electronic re- 



sponse function x a (uj) can either be computed from a micro- 
scopic model 43 or be determined from a fit to the Raman data 
using a phenomenological model function i 18 ' 40,4 1,45 These two 
approaches will be described in the following two subsections. 

It is instructive to note that formula (fTJ can be brought into 
the form of the widely used "standard" Fano profil e) 17 ' 18 ' 34 ' 46 



I P (uj) = C F (q + e) 2 /(l + e 2 ), 



with the asymmetry parameter 



q = -5(w)/[p(w)] 



(2) 



(3) 



and e = (uj — ujq)/T as in Eq. ([T). To extract renormalized 
phonon parameters the Raman spectra are often fitted using 
the above Fano profile, Eq. (O, with the Fano parameter q, the 
intrinsic phonon frequency ujo, and the renormalized phonon 
line width V (half-width at half-maximum, HWHM) kept fre- 
quency independent. While such an approach can give valu- 
able insights into the temperature dependence of the phonon 
lineshapes, the intrinsic electron-phonon coupling constants 
and the shape of the electronic continuum x(^) m the SC state 
cannot be determined. We will compare the simplified Fano 
formula (f2]i with our generalized theory [Eq. (Q~|i] in Sec.lIVI 



A. Phenomenological model of the electronic response function 

To apply the phenomenological model function [Eq. (Q]i] 
we need to either assume an expression that describes the real 
and imaginary parts of the electronic response function x( w )> 
or compute the electronic Raman response from a microscopic 
model. Let us first describe a phenomenological model func- 
tion for Following along the lines of Refs. 401 4711451471 
we express the imaginary part of the Raman efficiency with 
three terms 



p(uj) 



^w 2 + 



Ci 



Co 



1 



+ 4H 

(uj — > —Oj) 



— (u; — >■ — uj) 



(4) 



where e\(uj) = (uj — uji)/Ti and ea(o;) = (uj — wa)/^. 
The first term in Eq. (0]) models an incoherent electronic back- 
ground, which dominates the response in the normal state of 
cuprate superconductors. It is linear in uj at small frequencies 
and becomes constant for large Raman shifts 44 . The second 
and third terms are Lorentzians describing the pair-breaking 
peak located below 2 A, and the suppression of spectral weight 
at low frequencies, respectively. The latter is due to the open- 
ing of the SC gap. To reduce the number of free fitting pa- 
rameters, we set cl>2 = T2 = (uj\ — Yi)/2. The last two 
terms in Eq. (|4| decrease in intensity and shift spectral weight 
to lower frequency as the superconducting transition temper- 
ature T c is approached from below, {oj — > — uj) terms, similar 
to the Ci and C2 terms but with uj replaced by — uj, are es- 
sential to achieve the symmetry requirements for the Raman 
response. In the normal state, the electronic response func- 
tion is entirely described by the incoherent contribution [first 
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FIG. 2: (color online) B\ g Raman spectra of YBCOe.95 (A = 514.5 
nm) in the normal state at T = 100 K (a) and in the superconducting 
state at T = 20 K (b). Open circles show the experimental data, the 
solid curve the fitting result with formula Q} and four Lorentzians 
for the remaining A\ g phonons. The insets show the corresponding 
data after phonon subtraction. Panels (c) and (d) show the phonon- 
subtracted spectra for the xx and yy channels, respectively. The 
vertical lines indicate the maxima of the electronic peak intensity 
at 450 cm -1 and 480cm -1 , respectively. 



term in Eq. ©]. The real part of the electronic response func- 
tion, R(u>), is obtained from the Kramers-Kronig transform 
of p(uj). Since R(w) renormalizes the phonon frequency oj 
and the parameter So, formula (Q3 together with Eq. (0| and 
its Kramers-Kronig transform yield a self-consistent analysis 
of the Raman spectra. To compute the Kramers-Kronig trans- 
form of the incoherent part in Eq. (0}, a cut-off frequency ui c 
has to be introduced, which results in a constant offset in the 
real part of the electronic response R(oj). Provided oj c is cho- 
sen large enough, however, this error only leads to negligibly 
small corrections. 

We have employed a non-linear fit procedure with ten (six) 
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TABLE I: Extracted parameter values of the B\ g oxygen vibration 
in YBCOe.95 measured at 20 K and 100 K, respectively. Left part: 
asymmetry parameter So, intrinsic phonon linewidth (FWHM) 2To, 
and intrinsic phonon frequency ujo using the generalized Fano ap- 
proach of Eq.^T}. Right part: ujo extracted with Eq.© for different 
phenomenological electronic backgrounds [sqrt = first term in Eq.©; 
lin. = linear background with offset at uj — 0]. For comparison we 
display in the middle part loq obtained from inelastic neutron scatter- 
ing (INS) experiments by Reznik et al£^ 



independent fit parameters to the Raman spectra of YBCO6.95 
and YBC06.95:Ca in the superconducting (normal) state. C n 
and lot describe the intensity and position of the maximum of 
the normal-state electronic background given by the square 
root of a rational function. In the superconducting state 
C1/C2, Wi,2> an d Ti 2 describe amplitude, position and width 
of a Lorentzian function reflecting the region of the pair- 
breaking peak. Finally, r , lu$, and So effectively charac- 
terize amplitude, width, position, and asymmetry of a gener- 
alized Fano function (describing the B\ g phonon). Note that 
for not too strong frequency dependences of p(ui) and R(ui), 
the renormalization of the Fano formula due to the electronic 
background can be simply viewed as an offset of the param- 
eters entering in the Fano formula [cf. Eqs. (1), (2), and (3)]. 
The results of this fit procedure are shown in Fig.|2]and Fig. [3] 
together with the experimental data. Table U lists the corre- 
sponding fit parameters. Note that on general grounds the ex- 
tracted ujq and To parameter values need to be identical at a 
given temperature for all measurements of a given phonon. In 
Table U we show, however, the parameter values that corre- 
spond to the best (non-linear least squares) fit to the data. The 
differences reflect the error bars of our procedure. It is impor- 
tant to emphasize that the Fano profile of the £?i g -mode shown 
in Figs. [2] [3] and Ufa) results from the interaction with the 
electronic Raman signal; both electronic and phononic con- 
tributions and their interdependence can be described by Eqs. 
(Q} and (01). The final result agrees well with the measured 
data. Thus, our model allows us to some extent to disentan- 
gle the electronic and phononic parts of the Raman response, 
and to identify the shape of the electronic background. How- 
ever, strictly speaking, only a combination of parameters such 
as g 2 p and g 2 R can be extracted. In the next subsection this 
procedure will be improved for the SC state by employing a 
microscopic description of the pair-breaking excitations. 
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FIG. 3: (color online) Same as Fig.[2]but for the overdoped sample 
YBCOs.gsiCa. The vertical lines indicate the maxima of the elec- 
tronic peak intensity at 290 cm -1 (XX) and 355cm _1 (YY), respec- 
tively. 



B. Microscopic description of the superconducting state 

To describe the polarization-dependent electronic response 
function Xa(<^>) m the superconducting state, we employ a mi- 
croscopic model using a realistic tight-binding band structure 
with anisotropic hopping parameters and a superconducting 
gap with a mixture of d x 2_ y 2- and s-wave symmetries. Such a 
microscopic approach in conjunction with the analysis of the 
Raman spectra based on Eq. (Q~|i allows us to obtain precise 
information about the wavevector dependence of the super- 
conducting order parameter. In particular, we are interested in 
estimating the magnitude of a possible s-wave admixture to 
the d x 2_ y 2 pair wave function. 

The starting point of our calculation is the microscopic 
model introduced in Ref. |2lJ. The electronic Raman response 



function in the SC state for a given symmetry channel a is 
described by 16 - 21 . 43 . 48 



(5) 



9k(u>) being the Tsuneto function. The angular brackets de- 
note the average over the Brillouin zone, that is, 



V 



2T 



it) + 2E k 



irj — 1E\, 



(6) 



The quasiparticle dispersion E k = yej, 
fective one-band description of a single copper-oxygen plane 



A| contains an ef- 



£fc = -2t[(l + <y ) cos k x + {l-& Q )caak y ] 



-4f ' cos fc„ cos k„ 



and the superconducting gap 



(cos k„ — cos k, 



A.. 



(7) 



(8) 



For the computation of the electronic Raman response in the 
superconducting state, Eq. (O, we have assumed the same 
band structure parameters as in Ref. HI] together with A^ = 
30 meV and T e = 5 meV. 

The Raman vertices j a in Eq. (0 can be classified accord- 
ing to the irreducible representations of the symmetry group 
of the crystal. Since we shall consider a model with small 
distortions, S <C 1, A s <C A d , we use the notation corre- 
sponding to tetragonal symmetry, as in Sec. [TT] above. We can 
thus express, for example, the B lg Raman vertex as 



j Bi k oc t [(1 + 5 ) cos kx — (l- S ) cos ky] 



(9) 



To describe the combined electronic and phononic Raman 
response in the SC state, we use Eq. (Q]) together with the real 
and imaginary parts of the electronic Raman response in the 
SC state, Eq. ©. Furthermore, we assume that the bare fit 
parameters do not change as we go from the normal state to 
the SC state. 

In Figure Ua) we show numerical results obtained from 
our theoretical model and their comparison with the data on 
YBCO6.95 obtained in £>i g -polarization. The calculations 
have been performed at T = 20 K and for various s-wave 
contributions. The best description for both the £> lg -mode 
and the electronic response is found if A s is assumed to be 10 
percent of the maximum of the d x 2_ y 2-wave gap (short dotted 
line); this also accounts for the slight shift of the pair breaking 
peaks in xx and yy polarizations (Fig. [2] c and d). Further- 
more, our numerical results solely for the electronic response, 
i.e. the 2A-pair-breaking peak, are also displayed in Fig. [4] 
(dashed, dotted and dash-dotted lines). They are obtained by 
setting all phononic parts and the corresponding interference 
terms in Eq. (Q~|i to zero. Interestingly, we find that the pair- 
breaking peak shifts to lower energies with increasing A s , a 
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FIG. 4: (color online) (a) B\ g Raman spectra of YBCO6.95 in 
the superconducting state at T = 20 K (A=514.5 nm). Open 
circles show the experimental data. The short dashed, short dot- 
ted, and short dash-dotted curves are the theoretical result ob- 
tained with Eq. l[5}. The B\ g phonon was also taken into ac- 
count. We show results for various s-wave contributions A s /Ad = 
(dashed); 0.1 (dotted); 0.2 (dash-dotted). The dashed, dotted, and 
dash-dotted curves show the calculated imaginary part of the elec- 
tronic response, Xsi ( w )- (b) Subtracted spectra [20K (sc-state)- 
100K (n-state)] for the Bi g polarization channel of YBCO6.95. Be- 
fore subtraction the spectra were divided by the Bose factor. The 
solid curve depicts the theoretical result for 10% s-wave contribu- 
tion. 



fact which has been discussed in a previous paper by Schny- 
der et al.—. In addition, the cubic low-energy response, i.e. 
its (oj 1 2 A ) 3 -behavior found for A s = changes if A s ^ 0: 
we obtain linear correction terms which are, however, propor- 
tional to A s / Ad, and thus barely observable! 3 - 

Figure |4|b) shows the differences (x'$ — Xn)< which were 
obtained by subtracting the spectra at 20 and 100 K after di- 
viding by the Bose factor. The solid line is obtained after 
subtracting the model function (|4]i for the normal state [see 
Fig. |2 a)] from the results for the SC state when A s / Ad = 0.1 
[dotted line in Fig. Ufa)]. We find that the position of the 
pair-breaking peak at ~ 460 cm -1 (and partly its shape) is 
well described by our theory. This confirms that 10% s-wave 
contribution [short dotted line in Fig. Ufa)] yields an excellent 
description of the electronic Raman response. Furthermore, 
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FIG. 5: (color online) (a) B\ g Raman spectra of YBC06.9s:Ca 
in the superconducting state at T = 12 K (A=514.5 nm). Open 
circles show the experimental data. The short dashed, short dot- 
ted, and short dash-dotted curves are the theoretical result ob- 
tained with Eq. 1(5}. The B\ g phonon was also taken into ac- 
count. We show results for various s-wave contributions A s / Ad = 
0.05 (dashed); 0.15 (dotted); 0.25 (dash-dotted). The dashed, dot- 
ted, and dash-dotted curves show the calculated imaginary part 
of the electronic response, x'b 1 ( w )' (b) Subtracted spectra 
[12K (sc-state)-80K (n-state)] for the Bi 9 polarization channel of 
YBCOe.gsiCa. Before subtraction the spectra were divided by the 
Bose factor. The solid curve depicts the theoretical result for 15% 
s-wave contribution. 



the temperature-dependence of the related B\ g phonon is also 
reproduced (see below). Note that the difference between ex- 
perimental and calculated data at small energies is likely to 
be due to elastic impurity scattering which is not taken into 
account in our theoretical model4^ 

Finally, we return to the moderately overdoped sample. Ra- 
man data on YBCOg.95:Ca above and below T c in several 
polarization geometries are displayed in Fig. [3] In the ab- 
sence of specific information about the electronic band dis- 
persions of this material, we modeled these data by scaling 
the magnitude of the energy gap A^ by the ratio of transi- 
tion temperatures, keeping all other model parameters iden- 
tical to those used for YBCO6.95. The best fit was obtained 
for A s /Ad = 0.15 (however the estimated error bars are 
about ±0.05), slightly larger than the corresponding quantity 
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FIG. 6: (color online) Fano-analysis [Eq.l(2j] of three different 
phonons measured at a temperature of 20 K (A=5 14.5 nm). (a,b) The 
fit of the 340 cm -1 mode, which corresponds to the out-of-phase vi- 
brations of the planar oxygen, along yy and xx polarizations, respec- 
tively.(c,d) same as (a,b) but for the in-phase planar oxygen mode 
around 440 cm -1 and apical oxygen mode at 501 cm -1 . Open cir- 
cles show the experimental data. The dash-dotted lines are the results 
of the overall fit to the experimental data. The solid lines are the re- 
sults of fits to Fano-profiles [Eq.©] as described in the text. The 
dotted lines correspond to a linear background. The intensity units 
are arbitrary but the same in the four vignettes. 



in YBCO6.95. Note that the quality of the fit is comparable to 
the one for YBCO6.95, although separate plane and chain sub- 
systems were not introduced in the analysis^ The previously 
observed difference between spectra in xx and yy geometry 
was confirmed (Fig. [3]), but our model calculations suggest 
that this is a consequence of the s-wave admixture to the gap= 
obviating the need to introduce quantum interference between 
scattering from chains and planes^2i£i 



IV. TEMPERATURE DEPENDENT PHONON 
LINESHAPES 

A. Anisotropic Fano profile 

Several phonons in YBCO display a pronounced asymmet- 
ric lineshape suggesting a strong interaction with the elec- 
tronic continuum. As shown in Fig. |6]the asymmetry is most 
pronounced for the 340 cm -1 mode. Interestingly, the phonon 
peak reflecting the vibration of the apical oxygen at 50 lcm -1 
exhibits a strong asymmetry for a polarization of incident and 
scattered light along the a-axis, whereas along the fe-axis the 
phonon appears to be almost symmetric (see Fig. |6j. Two ad- 
ditional modes are present in this spectral range. The mode 
at 440 cm -1 originates from an in-phase vibration of the oxy- 
gen atoms 02 and 03. Additional modes are present at about 
472 cm" 1 and 480 cm" 1 for polarizations parallel to the a- 
and 6-axis, respectively. As mentioned above, these modes 



are Raman forbidden, but correspond to IR allowed vibrations 
involving the Cul-Ol chains^ Due to defects (oxygen vacan- 
cies) in the Cul-Ol chains they become Raman active. 

Similar to previous temperature dependent Raman experi- 
ments on YBC0 7 _,?^ and HgBa 2 Ca 3 Cu40io+,5ri we have 
fitted the phonons of YBCO6.95 by using simple Fano pro- 
files. The solid lines in Fig. [6] are the results of fits to the 
experimental data using Eq. (O , the dash dotted lines corre- 
spond to the resulting fitted lineshape of the entire spectrum. 
The calculated profiles agree well with the measured spec- 
tra. In order to obtain estimates of the intrinsic phonon posi- 
tions and linewidths, we have corrected the peak position and 
linewidth by the procedure used in Ref. [2CJ We find for the 
xx-symmetry, i.e., polarizations along the a-axis: 



w 



r/ gee 



and 



-pa; a: 

1 FWHM 



2r 



(1 



g 2 



(! -llx) 



(10) 



^fwhm denotes the intrinsic full width at half maximum 
(FWHM) 5 ^. For the yy-polarizations one obtains an analo- 
gous equation. 

Fits to the simplified Fano profile of Eq. (O are much less 
complicated than the fitting procedure to the full spectrum dis- 
cussed above. This is a key advantage especially in situations 
in which several closely spaced phonons partially overlap as is 
the case, for instance, for the apical-oxygen vibration in Fig. 
[6] The multi-parameter, global fit [Eq. (Q~|)] yields unstable re- 
sults for these phonons. It is important to note, however, that 
the quantities extracted from simple Fano fits are renormal- 
ized by the electronic response function, and therefore deviate 
slightly from the intrinsic phonon frequency uj q and the bare 
linewidth Tq of Table Q] We will compare the results of both 
procedures in detail at the end of the next subsection. 



B. Superconductivity-induced changes in the position and 
linewidth 



Figure [7] shows the temperature dependence of the ener- 
gies and linewidths of two particular phonons, the 340 cm -1 
and 501 cm" 1 modes, measured with light polarization along 
the crystallographic a- and fe-axes. The spectra were taken 
at temperatures ranging from 20 to 300 K. The tempera- 
ture dependence of the phonon energy and linewidth in the 
normal state arises from anharmonic phonon-phonon interac- 
tions, i. e., the decay of a high-energy optical phonon into two 
phonons of lower energy with opposite momenta For 
simplicity, assuming the resulting phonons to have the same 
energy 56,57 and using Bose-Einstein statistics, this decay pro- 
cess leads to r onfc .(T) = rj= h °(l + 2n(w p /2)). This pro- 
cess, implying decay through real transitions, is strictly valid 
for the line width, but has been also used for the frequency 
shift although, in this case, virtual transitions also play a role. 
In Fig. [7] both linewidth and peak position were fitted si- 
multaneously for the temperature range above T c (solid lines). 
Both of these quantities show abrupt changes at the SC tran- 
sition temperature due to the opening of the superconduct- 
ing gap, as previously observed in YBCO (see for example 
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FIG. 7: (color online) Frequency and line width (FWHM) ob- 
tained using Eqs. ((2) and dlOt versus temperature of the 340 cm -1 
mode and the 501 cm -1 mode of YBCO6.95 for the xx- and yy- 
polarization, respectively (A=514.5 nm). Below T c a comparison 
between calculations for s — A 3 / = (dashed line), 0.1 (dash- 
dotted line), and 0.2 (dotted line) is displayed. The solid line cor- 
responds to the temperature dependence of pure phonon-phonon in- 
teraction (interpreted as Klemens decay into two phonons of equal 
frequencies). 



Refs. [33lf53l) and also other superconducting compounds like 
HgBa 2 Ca 3 Cu 4 Oio^ and Bi 2 Sr 2 CaCu 2 8 ^ i S 

Within the error bars of ~ 0.8 cm" 1 in (a)-(c) and 
2.0 cm -1 in (d), the phonon peak position and linewidth of 
1 mode is the same for the xx and yy polar- 



the 340 cm" 
izations (see Fig. 



|2h and b). Taking the phonon positions 



and linewidths in the SC state as uj p and T s p , the maximum 
change in the phonon position and linewidth is obtained from 
Acj = oj s p {T)- uj p {T = 100 K) and Ar = T s p (f)- 
T p (T = 100 K)r&. where T denotes the temperature at which 
the maximum of the SC-induced changes occurs. The addi- 
tional softening below T c due to the electron-phonon interac- 
tion is about -6 cm -1 . The change in linewidth (FWHM) is 
+3 cm -1 , reflecting a broadening. For the vibration of the api- 
cal oxygen, i.e. the 501 cm -1 mode, the corresponding shift 
and broadening are +3 cm -1 and +5 cm -1 , respectively. The 
SC-induced changes for the 340 cm -1 and 501 cm -1 modes 

are in good agreement with previous data obtained on twinned 
YBCQ^eo 

The SC-induced changes in the phonon linewidth and peak 
position (for a phonon labeled as fi) can be related to changes 



r) — > 0), resulting from the interaction between the phonons 
and the electronic system below T e i 21 - 48 ' 59 The induced fre- 
quency shifts Auj are related to the real part of the phonon 
polarization II by 48 ' 61 ' 62 



Aw 

U) 



1 



AT(0) 



A Ren, 



(ID 



while the induced changes in the linewidth can be calculated 
via 



AT 
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(12) 



Here, N(Q) denotes the electronic density of states at the 
Fermi level and A = 2 £ k £ p / ^ | 5 £ | 2 F k »<5(e k ) 
is the dimensionless electron-phonon coupling constant. 
F£(uj) denotes the spectral function for phonon \i un- 
der consideration.*- Taking into account screening effects, 
i.e. the long-range Coulomb force, the self-energy is given 

jj., . 21,48,61,62 



// \2 \ \.9fc(AM' 

E» = -((<„) W) + ^ (W) ' . 

where the angular brackets are defined by Eq. (O. The sym- 
metry of the optical phonons is reflected in the matrix element 
g£ Q . The electron-phonon coupling of phonons of A lg and 
B lg symmetry are in a first approximation given by 



9kfi = 9 Blg ( cos K ~ c osk y )/2, 
9k!o = 9a 1cj (cos k x + cos k y )/2. 



(14) 



in the phonon self-energy E (u>) 



n(w + ii]) (with 



with the electron-phonon coupling strength g B and g A 
(Note that there is more than one phonon of A\ g symmetry). 
In general, the phonons below the energy of the superconduct- 
ing gap 2A mQX should soften below T c (i.e. they should shift 
to lower energies), whereas phonons above 2A maa; should 
harden. This is confirmed in our experiments (not all data 
are shown here) and in previous work (see Refs. [l9l [33L and 
l53l) . The energy of the apical oxygen phonon (501 cm -1 = 
62.5 meV) is right at the gap energy and is therefore sensi- 
tive to small changes in the energy of the SC gap. The 340 
cm -1 mode is well below the SC energy gap for T >C T c . 
With increasing temperature the energy of the 2A gap shifts 
to lower energies and moves through the energy of the 340 
cm -1 mode. This explains the maximum of the linewidth at 
75 K. A similar behavior was observed by Limonov et alM. 

Finally, we discuss the role of the s-wave contribution 
to the superconducting gap which has been introduced in 
Eq. ([8). Our numerical results obtained with Eqs. ©-(O and 
Eqs. dTTb— (TT4T> using the fits reported in Sec. [Ill] are also dis- 
played in Fig. [7] We compare A s = (dashed line) with 
A s = 3 meV (dash-dotted line) and A s = 6 meV (dotted 
line). We find that the results obtained with 20 percent s-wave 
contribution cannot describe our data. This is clearly visible 
in Figs. |3 a) and (c) which show the SC-induced changes Alu 
in the position of the corresponding phonon. In particular, for 



9 



the 501 cm _1 -mode which is known to be very sensitive to the 
superconducting gap, 2 - one predicts a softening for A s = 6 
meV, while instead a hardening is observed. Therefore, our 
data imply an upper limit of A s / Ad = 0.2. The best agree- 
ment is obtained for 10% s-wave admixture. 

At the end of this subsection, we are contrasting our gen- 
eralized Fano theory [see Eq. ((TJ] with the standard Fano ap- 
proach described by Eqs. © and ( [Tol l. The main difference 
between both approaches is the theoretical description of the 
electronic Raman response. While in the standard Fano ap- 
proach the background is assumed either to be linear (see 
Fig. |6]l or to follow a square-root behavior— [first term in 
Eq. in both the normal and SC-state, our generalized Fano 
theory is able to take the rearrangement of spectral weight 
due to the opening of the superconducting gap into account. 
Within our microscopic description [Eqs. ©-(HJl] it is then 
possible to determine the ratio A s /Ad- Another difference 
between both Fano theories concerns the asymmetry parame- 
ter q: in our generalized theory it becomes w-dependent which 
allows a self-consistent description of the £> l9 -phonon and 
the electronic Raman response from 80 to about 1000 cm -1 . 
Thus, in other words, no specific frequency interval close 
to the phonon peak position [see Fano lineshape analysis in 
Figs-ISta) and (b)] needs to be selected. 

We have summarized our comparison in Table I in which 
the intrinsic frequency of the B\ g -phonon extracted from both 
Fano theories is displayed: ujq obtained from Eq. ((TJ [left part, 
third column], uiq obtained from Eqs. (O and ( TTOb either for 
a square-root [wo( s qrt)] or linear background [wo(lin.)] (right 
part). These values have to be compared with uio obtained 
from inelastic neutron scattering experiments [ujq (INS)] (mid- 
dle column). One clearly sees that the values for uiq ob- 
tained from our generalized Fano theory are similar to those 
measured in INS experiments. On the other hand, the val- 
ues extracted from the standard Fano approach differ by 1.1 
[Bi 3 (100K), cj (sqrt)] to 5.9 [XX(20K), w (lin.)] wavenum- 
bers. Note, however, that these values can be substantially 
modified if a specific frequency interval around the peak posi- 
tion of the Bi 9 -phonon is defined and the corresponding line- 
shape analysis is then restricted to this interval (see Fig. [6]) 
These improved values are shown in Fig. Q In the case of 
the apical oxygen vibration analyzed in Figs. |6jc) and (d), a 
fit with Eq. ([T]i is hardly possible due to the multiple peak 
structure between 400 cm -1 and 540 cm -1 . Therefore, we 
have used Eqs. (|2), (0), and ( [Tol l. On the whole, we find char- 
acteristic differences between our self-consistent, generalized 
Fano theory and the standard Fano approach, but our conclu- 
sions about the ratio A s /Ad are independent of the theory 
used provided that Eqs. (|2), (0), and ( fTOb are restricted (or 
strongly weighted) to an interval around the peak position of 
the phonon being investigated. 



C. Temperature dependence of the asymmetry parameter 

Figure [8] shows the temperature dependence of the asym- 
metry parameters and q~y in Eq. (O as obtained from the 
fits of the Fano profile to the experimental spectra. We discuss 
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FIG. 8: (color online) Fano asymmetry parameter (A = 514.5 nm) 
of both polarization channels for the 340 cm -1 and 503 cm -1 mode 
of YBCO6.95, respectively. Results of the theoretical calculations for 
different s-wave admixtures (0, 0.1, and 0.2) are also shown. 



these data in the framework of the approximate expression 

_ To Ha 



9<?Xo 



(15) 



tors . J 'i-"' t " We first focus on temperatures above T c and note 
-0.12 for the 340 cm -1 mode, whereas 



already applied in Raman work on other high T c superconduc- 

37.54.65 

that q~l « q-y 1 
a pronounced difference in the asymmetry parameters along 
in xx and yy geometries is observed for the 501 cm -1 mode 
(ixx ~ —0-20, q~y ~ —0.08). As g a is a materials parameter 
independent of the light-field configuration and is weakly 
cj-dependent, the strong variation of asymmetry parameters 
must be attributed to differences in T a /j a . Our results are 
in fair agreement with the phononic and electronic Raman ef- 
ficiencies (proportional to \T a \ 2 and \~fo\ 2 , respectively) cal- 
culated in refs. |66|67| using density functional methods. For 

instance, the calculated l 7 ^!"" (~ ^-) in the normal state is 

about 1.2 and 0.65 for the 340 cm -1 and 501 cm^ 1 modes, 
respectively. The changes of the asymmetry parameters with 
temperature in the normal state may reflect a temperature de- 
pendence of 7c,, i.e. parameters in the electronic band struc- 
ture as well as that of the anharmonic linewidth, but their de- 
tailed origin cannot be disentangled 

Moving to the superconducting state, we note that the only 
parameter in Eq. (TT3T > expected to change significantly across 
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T c is the response function \"- The temperature depen- 
dence of q~]; and q~y below T c therefore also reflects the ab- 
anisotropy in the 2A-gap. Thus we have computed the tem- 
perature evolution of the asymmetry parameters in the super- 
conducting state, using the same model parameters described 
in SecHIIl The reasonable agreement with the data (Fig. [HJ 
demonstrates the self-consistency of our approach. However, 
since the superconductivity-induced modification of and 
q~y is relatively subtle (compared to that observed in the nor- 
mal state), this analysis does not provide constraints on the 
ratio A s /A,i beyond those already discussed above. 

It is interesting to note that other A g modes also reveal a 
temperature-dependent asymmetry. In particular, we have ob- 
served that the Ba mode near 115 cm -1 is symmetric {i.e. 
q~ x = 0) for T — > and its asymmetry increases mono- 
tonically to q^ 1 = 0.25 (not shown). Its intensity, however, 
is rather weak compared with the B\ g mode (and this applies 
even more so to the Cu mode near 150 cm -1 ), so we have not 
used these modes further in our discussion. Their asymmetry 
has been analyzed in more detail in Ref. 69, however without 
considering a possible s-component to the <4;2_„2-wave gap. 

V. SUMMARY AND CONCLUSIONS 

In this work we have revisited the superconductivity- 
induced effects on phononic and electronic Raman scatter- 
ing in detwinned slightly overdoped YB^CuaOg.gs and mod- 
erately overdoped Y0.85Ca0.15Ba2Cu3O6.95 single crystals, 
both from an experimental and a theoretical point of view. In 
particular, we have performed a detailed study of the in-plane 
anisotropics in the electronic continuum and in the phonon 



lineshapes and assessed their implications for an s-wave ad- 
mixture to the d-wave superconducting gap. To this end, we 
developed and applied a formalism that treats the frequency 
dependence of both electronic and phononic Raman scatter- 
ing on equal footing. Thus we can disentangle both parts 
and clarify the role of interference terms. Due to this pro- 
cedure we are able to extract the intrinsic frequency ojq of the 
phonon position reflecting the intrinsic electron-phonon inter- 
action. A comparison shows that ujq agrees well with INS data 
at the T-point. The best agreement with the Raman data was 
obtained by model calculations based on admixtures of 10% 
and 15% ± 5% s-wave contribution for YBa2Cu.3O6.g5 and 
Y0.85Ca0.15Ba2Cu3O6.95, respectively. This agrees with val- 
ues obtained by other experimental method a 3 i 4 ' 5 ' 6 ' 7 i 8 and con- 
firms the previously observed trench of an increase in the s- 
wave contribution with increasing doping level. Our data do 
not show evidence of the previously reported^ unusual quan- 
tum interference between electronic Raman scattering from 
planes and chains in overdoped Yi^Ca-jT^CusOy-^. 
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